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Detecting compact dark matter in galaxy clusters via 
gravitational microlensing: A2218 & A370 
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ABSTRACT 

After decades of searching, the true nature of dark matter still eludes us. One poten- 
tial probe of the form of dark matter in galaxy clusters is to search for microlensing 
variability in the giant arcs and arclets. In this paper, a simple method is introduced 
to characterize pixel variability in the limit of high optical depth to microlensing. 
Expanding on earlier work, the expected microlensing signal for two massive clus- 
ters, A2218 & A370 is calculated. It is found that the microlensing signal depends 
sensitively upon the mix of smooth and compact dark matter in the cluster. Com- 
parison of two deep exposures taken with James Webb Space Telescope or two hour 
long exposures taken with a 30-metre class telescope in two epochs separated by a 
few years will possibly detect about a few dozen pixels which show strong variability 
due to microlensing at 5er level, revealing wealth of information on the microlensing 
population. 

Key words: gravitational lensing : microlensing : high optical depth - dark matter 
: compact - galaxies : clusters : individual (A370, A2218) 



1 INTRODUCTION 

There is significant evidence that most of the energy den- 
sity in the Universe resides in the forms yet unknown, but 
the physical nature of this dark matter is an issue of live 
debate and is far from being understood lILahav fc L iddlc 
l2003l :lMunos 2003). Typical explainations invoke either com- 
pact objects, such as primordial black holes, stellar remnants 
or planets, or continuous material such as weakly interact- 
ing massive particles, unseen hydrogen or more exotic ex- 
plainations. Potential clues to the nature of dark matter 
have recently been uncovered with the detection of neutrino 
mass and observations of compact massive objects thought 
to be mostly white dwarfs in Galactic neighbourhood (e.g., 
iFukuda et all2000llAl"cock et all200ch . However, given their 
mass neutrinos cannot contain mor e than ~13 per cent of the 
dark matter mass budget (iHu. Eisenstein fc Tegmark 1998). 
Furthermore, the halo microlensing programs constrain com- 
pact objects in the mass range of (10~ 7 — lO -1 )./ ^ to mak- 
ing n o more than 10 per cent of the Galactic halo JSadouletl 
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Gravitational lensing has allowed the detailed recon- 
struction of the projected ma ss distribution in galaxy clus- 
ters fe.g.. lEbbels et al]ll9 98'l. Unfortunately such analysis 
does not probe the fundamental nature of dark matter. How- 
ever, although deflection angles due to compact objects in 
the dark matter are negligible their derivatives may be sub- 
stantial which means that although the granularity of dark 
matter cannot change the morphology of lensed systems, this 
microl ensing can have significant i mpact on the observed 
fluxes (Liebcs 1964; Paczvriski 1986). As the mutual config- 
uration of sources and lenses change with time, this results 
in variability of the observed flux. 

Several microlensing programs have been proposed 
and undertaken, the most successful of them being 
the microlensing experiments towards Galactic bulge 
and Magellanic Clouds by MACHO EROS and OGLE 
collaborations l^ coc^^a^ 200C ; lLasserre et al.l .2000; 
lUdalski. Kubiak fc Szvmanskil 1997 1. Another approach is 
pixel lensing when the search focuses not on individual 
sources but on a large number of stars seen as a single pixel of 
the image. Events of strong lensing of a single star can be de- 
tected in the pixel light though they never dominate it. Such 
events have been discussed in c onnection with the galax- 
ies M 31 and M87 (ICrott dll99sl : iBaillon et alJll993tlGouIdl 
1995), and a few observational programs monitoring these 
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galaxies have been implemented iCrotts fc Toma nevlll996t 
lAnsari et aHll997tlRiffeser et alJl200ltlBaltz et al.ll2003t) . 

Further away from the Galaxy, it has been pro- 
posed to search for microlensing-indu ced variability of 
quasars seen through galaxy cluster s JWalker fc Ireland! 
Il995t iTadros. Warren fc Hewetdll998h . This has been un- 
dertaken for the Virgo cluster, although its close proximity 
ensures that the optical depth to microlensing is low. Al- 
though analysis of the quasar vari ability based on a long ob- 
servation series bv lHawkinj <ll996f) let him conclude that the 
dark matter may be dominated by Jupiter-mass microlenses, 
this idea is not widely supported. One of the difficulties here 
is that quasars are exp ected to possess intrinsic variability 
ijZackrisson et al.ll2003TI which significantly complicates mi- 
crolensing studies. Totani (2003) proposed to explore the 
advantage of the recent discovery of a galaxy cluster found 
just behind the rich cluster Abell 2152 but this has yet to 
result in an observational program. 

Another possibility is to look for the variations in sur- 
face brightness of strongly lensed distant galaxies - e.g., gi- 
ant gravitationally lensed arcs (Lewis & Ibata 2001; Lewis, 
Ibata & Wyithe 2000). They are found in the regions of 
large magnification and therefore high optical depth. As- 
suming that microlenses make up an appreciable fraction of 
the lensing mass, it means that in any instant all the stars 
in a pixel are subject to strong microlensing and therefore 
an investigation of the surface brightness variability is ef- 
fectively reduced to an investigation of the behaviour of the 
sum of a large number of individual fluxes with fairly well 
known individual statistical properties. ^^^_^^^^ ^_^ 

This study extends the analysis of iLewis et al] (2000) 
and provides a general method for the analytic calculation 
of variability patterns given the convergence and shear val- 
ues. We apply it to the two well-studied galaxy clusters - 
Abell 370 and Abell 2218. In the next section we describe 
and justify the method we use, and then calculate the value 
of individual microlensing-induced variability dispersion as 
a function of standard microlensing parameters. Section 3 
discusses observational prospects for detecting this sort of 
variability. We discuss our results in section 4. 



time scales greater than typical individual variability scale 
given by the time of Einstein- Chwolson radius crossing 
(r e /v), and therefore are interested in the probability distri- 
bution function of the summed pixel flux. The typical fluxes 
in pixels forming images of strongly lensed galaxies corre- 
spond to the range in luminosity of thousands to billions of 
solar luminosities - and therefore contain very large numbers 
of stars. This fact along with the variability independence 
immediately suggests using the central limit theorem to infer 
the statistical properties of the pixel flux. 

Indeed, let the intrinsic (unlensed) fluxes of the pixel 
population stars be {Li},i = L.A^, N being the number 
of stars in the pixel. If we neglect the fraction of stars in- 
trinsically variable at the same level and time scale as the 
variability caused by the microlensing, the role of Li is to 
normalize the variability in magnification factor fa arising 
as the source stars move through the magnification map: 



Lobs r 
i — J^ifa 



(1) 



and fit can be considered a random variable. 

The probability distribution of \ii does depend on the 
individual characteristics of the source - mostly, the size of 
its disc. However, as we will see in the following section, 
this dependence is not strong and, to show the validity of 
Gaussian approximations it is enough to assume that fa dis- 
tribution is identical for all i and only depends on (macro) 
lensing parameters in the pixel. 1 

The flux observed in the pixel is given by 



Lobs I 
\{Li} 



N 

E 



Lifii 



and its average value is 



L obs 



{l,} =M 



N 



Li, 



(2) 



(3) 



where the bar denotes the averaging over the lensing config- 
uration at a given convergence and shear (assuming a sort of 
ergodic hypothesis this is equivalent to averaging over time). 
Let us define 8 Li = (fa — 71) Li and consider the devia- 
tion of L obs \{ L .y from its average value: 



2 GENERAL METHOD 

2.1 Statistical approach 

The goal of this section is to learn to characterize the vari- 
ability in pixels of images of distant galaxies with a few 
parameters. We show that the high optical depth and large 
numbers of stars forming the pixel allow one to reduce this 
to a single parameter - dispersion - as the observed flux 
distribution can be approximated by the Gaussian function 
with good accuracy. 

The major assumption we make in our investigation is 
that the fluxes of individual stars, which vary as stars and 
lenses move with respect to each other, do so independently. 
This is a natural assumption since individual stars in the 
source plane are distributed in a random manner. It might 
be worth considering the presence of some form of agglom- 
erations in pixels but this seems to be an unnecessary com- 
plication at the present stage. 

We consider statistical properties 'as static' - i.e. on 



sl \{l,} ='y]Li(pi-n) 



(4) 



Let us also define the second and third moments of the dis- 
tribution in fi: aZ, = (jJ> — 7i) 2 and = \fi — 7i| 3 . Clearly, the 
corresponding moments for 6 Li equal Lfa^ and Li/3. 

According to Lyapunov theorem, the actual probability 



distribution function F of SL\ 
approximation 

1 _«i 2 



{Li} 



tends to the Gaussian 



$(5L) = 

x/2^2 

where we dropped {Li} subscript for clarity, with 



°l £ L * 



(•») 



(6) 



We also assume that these parameters are constant throughout 
the pixel and this is a natural assumptions in most cases. 
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The accuracy of this approximation in Kolmogorov (£00) 
measure p(F, "J>) = sup \F(x) — $(a;)| is not worse than cA, 

x£lZ 

where 



A 



^(E^?) 3/2 



(7) 



x/lO+3 
6V2tt 



0.41 



jBerrvl flQ4ll; lEsseenl Il942l) and c 
l|Chistvakovll200lft . 

Although the value of the second fraction in A may, 
depending on {Li}, be as large as unity its typical value 
is of order (L 3 )/((L 2 } 3 ' 2 V7V), where the angled brackets 
denote averaging over the luminosity function (which, up to 
a distance- and band-dependent constant, is the distribution 
of intrinsic fluxes) . 

With the lum inosity function given by 
IJahreiss fc WielenUl997l) 2 . (L?) « 10 3 L| and (L 2 ) w 10L; 
therefore 



0> 



cA^ 10 



P 1 

" 3 VA 7 ' 



(8) 



To estimate the values of /3 and a^ it is sufficient to use a 
rather coarse 'model' probability distribution density p((i) 
which is normalized to unity and has three basic properties 
established theoretically: 



(i) p(fi) = at n < 1 

(ii) J dfip(p)p - Hth, IHh - 1(1 - 

(iii) p(p) ~ p~ 3 at ^ » 1. 



k J -7 



The job is done by the following 'model' p(/J.)'- 
2(p -lf(p-l) 



p(m) 



[(m-i) 2 + (mo-i) 2 ] 2 



(9) 



The second condition implies fi — 1 = 2{fj, t h — 1)/tt. 

This distribution does not possess the second moment, 
let alone the third, as a result of the (iii) property. How- 
ever, the finite size of the source places a cut-off p ma x at 
the high values of /x (for a single point mass jj, m ax is nearly 
inversely proportional to the source size as was shown by 
Liebcs (1964)). Since fj, ma x ^> Po it does not affect ei- 
ther norm or the first moment and therefore with <r 2 w 

2plln(fJ, ma x/po) ~ Pt h ln(TTHmax/2pth) and (3 « Hth^max 

we have 



cA; 



10 



JV M*h l n (^fJ-max/2p t h) 



(10) 



With typical for gravitationally le nsed arcs values of pith ~ 
10, ^ maa; ~ 100 (Lewis et alJl2000|) and N ~ 10 4 - 10 6 in a 
pixel cA < 0.01 — 0.1. This means that when talking about 
a deviation of at least one standard value a, for which the 
Gaussian probability is w 0.15, one can be sure that the 
actual probability of such a deviation is not less than (5 - 
25) per cent. 

Strictly speaking, the minimum magnification value 
for microlensing at high optical depth is greater than 
one used in (i) as was shown by ISchneid er (1984), and 



2 The work of Ijahreiss &: WielerJ J1997TI presents V-band lumi- 
nosities which are of rather limited interest for exact cosmologi- 
cal predictions where K-corrections are to be taken into account; 
however, these numbers provide sensible estimates for the quan- 
tities under consideration. 



one could rather use some m odel value for this quantity 
iJBartelman fc Schneiderll99Cl) . However, this does not have 
much impact on the estimate of the validity of our approxi- 
mation. Perhaps more important is that due to the value of 
minimum magnification which is greater than unity, Gaus- 
sian approximation clearly cannot hold exactly as it as- 
signs non-zero probability to flux values below the minimum. 
However, this inconsistency is well inside the uncertainty of 
our method, given by 1101 and does not affect our results. 

The initial task is thus reduced to calculating the only 
parameter of a centered Gaussian distribution - its disper- 
sion. As {Li} is not known a priori (and neither it can be 
known well a posteriori), 8L\{L i y is to be averaged over all 
possible {Li}. This can easily be done by considering the 
following three random variables: 

AL = AL° + 5L (11) 

where AL = ^G^-J^I*}), AL° = /J^( L »^( L »)) and SL 

i i 

is the value of interest. Here, again the bar denotes averaging 
in the /j, domain while angled brackets mean averaging over 

{Li}- 

As these three quantities are (nearly) Gaussian and un- 
corrected (the correlation vanishes when averaging over /x), 
the following relation holds: 

"SL = "AL - °~AL° ■ (12) 

Clearly, 

"L=((^(i-(i.))) 2 )=^L (13) 

where g 2 is the dispersion of individual flux 
jJahreiss fc Wielenlll997li . 

For the second quantity we may write: 



2 

°~AL 






(pL- p{Li)) 



(14) 



= N(p 2 L2~2pLiL l {Li)+JI 2 {L i y) 

= N^LD-jfiLi) 2 ) 

= N(o-l(L 2 t )+Jl 2 a 2 L ) 

and therefore 

ah = N{LJ)al . (15) 

The dispersion of the quantity SL/(N~p{Li)) - the relative 
fluctuation - is thus given by: 



2 _ 1 (Lj) o- 2 



N {Li) 



(16) 



As iV is not known either we will just divide the observed 
flux of the pixel L t, s by the mean magnification factor and 
mean individual stellar flux to get a first-order estimate: 

Lobs 



N = 



p{L t ) 



Thus, 



SSL = 



{Lj) 

{Li)L b 



'sin, 



(17) 



(18) 



with the first factor in this formula being approxi- 
mately 6.02A/L o i, s /L(T) for the luminosity function of 
Ijahreiss fc WielerJ (I1997T) and the variability extent 
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F 2 = Zjl - E. _ i 

L H — —2 —2 ' 

p IX 

This quantity is therefore of our prime interest. 



(19) 



2.2 Evaluation of e 2 

In calculating the var iability e xt ent e 2 we em- 
ploy the method of INeindora J2003fl . who im- 
proved and ge neraliz ed previous wo rks o f 
iDeeuchi fc Watsonl Jl987h. ISeitz fc Schneider] (Il994l) 
and ISeitz. Wambseanss fc SchneideJ ill 994f) to make possi- 
ble the calculation of microlensing correlation functions in 
the case of non-zero shear. We, however, slightly modify his 
equations and evaluation method for our specific needs. 

Let z and C, be the light ray positions in the lens C and 
source S planes, respectively. The normalized lens equation 
is then (Kayser, Refsdal, Stabell 1986; Paczyriski 1986) 



C = JoZ — sign(l — « c ) \ ? 

i 

where 



Z — Zi 

Iz-zip 



Jo 



1+7 






1-7 



(20) 



(21) 



Here 7 = 7'/(l — k c ), while n c is the smooth matter cover- 
gence and 7' is the shear, both expressed in critical units 



4ttGD 



E = 
and 

r =_^_ 

" 4GD ' 
where D is a reduced (angular diameter) distance 

DlsDql 



D = 



Dos 



(22) 



(23) 



(24) 



The masses of the microlenses rrij are given in units of M 

- the quantity which also defines Einstein radii z and C,o - 
physical length units in which z and C of 12011 are expressed 

- in lens and source planes: 



4GMo 



DlsDql 



1 — K c 



Dos 



/4GM 0|1 ^DlsDos 

Co = \l 11 — K c 



Dol 



(25) 



(26) 



In the case we consider the microlensing shear S(z) = 
sign(l — k c ) ^mi(z — Zi)/|z — Zi| 2 - the second term in 120H 
- is an isotropic random variable. Changing the sign of 7 
has the effect of only redefining coordinate axes and since 
this is not of interest for us we drop the sign(l — k c ) factor 
in I20H and use absolute value of 7 from now on. 

The magnification factor at a p oint C in the sou rce plane 
may be written in an elegant form (Ncindorf 2003): 



KC) 



1 — K c 



S(C~ J D z + S(z))d 2 z. 



The average of jj, is 

1 
Hth 



|[1- Kc ] 2 [(1-k) 2 -7 2 ]| 



(27) 



(28) 



where n = 7r(m)n is the scaled microlensing optical depth 
with n being the surface number density of microlenses and 
angled brackets now meaning averaging over microlens mass 
distribution (j>(m). 

We consider a Lambert disc - a disc of uniform surface 
brightness - with radius R in ( units and total flux I . The 
average value of observed flux (J) does not depend on the 
microlens mass distribution and is only a function of lensing 
macro parameters: 



(I) = Io^th 



with 



J = Jo — «1 



h 



(1-K c ) 2 det J\ 



l-K+7 

l-«-7 



(29) 



As shown in the Appendix A, the following relation holds 
for e 2 : 



2 _ (J 2 ) _ 2|detJ| 

£ " + (7) 2 7T.R 2 



ds'drdx— x (30) 



'+t44" 



[na(s ,x)~t-* s cosx] r /2 



The integration is done from zero to infinity in p and s' and 
from zero to 2-k in ang ular variable y. Function a(s, \) ln 
this formula is given bv INeindora (|2003j) and is described in 
the Appendix A. 

Introducing the function 



B(s, X ) 



na(s, x) + is cos x 



,s 
a{s,x) 



(31) 



Tv{m) 



+ i cos x 



and performing the variable change r — > p = r 2 s/2, $' 
s = R 2 s' /2 we get the following integral to be evaluated: 



a _ jdetJ| 



dss x 



(32) 



d P j 2 (Vpi) Jo ( 7 p) /d xe - ps (iM 



This evaluation is done in Appendix B under the assump- 
tions that R < 10~ 2 - 10~ 3 < 1 and 7? < 7. 

The former is plausible, since in a cosmological situation 
with Z) ~ 10 2S cm, the length scale I26H is 



Co ~ 10 1 



M 

Mr, 



1/2 



1 1 1/2 

1 — k c \ ' cm, 



while typical physical sizes of sources are ~ 10 



10 1 



However, it immediately places a constraint on the microlens 
masses and the smooth component convergence: 



M \1-k c \ ^1O" 4 M 



(33) 



Therefore the results derived below are not directly applica- 
ble to situations where 1331 is not fulfilled, which may be of 
interest when Jupiter-mass lenses are involved or for detailed 
investigations of microlensing in the region \k c — 1| ^C 1. Mi- 
crolensing of large sources was considered numerically by 
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Figure 1. The contour lines of e M as a function of convergence 
in microlenses re and shear 7. 



iRefsdal fc StabelH (Il99lf) in the ze ro shear case and later in - 
cluding the effect of a shear term iRefsdal fc Stabelllll997t) . 
Quite naturally, increasing the size of the source suppresses 
microlensing-induced fluctuations, averaging them over less 
correlated regions of the magnification map. However, for 
sources as large as R ~ 30 they find values of e M ~ 0.1 in 
a range of $J K ^ 2 and ^ 7 Sj 0.4 (it was found that 
e^ w 2k/ R 2 in zero shear case). This is only an order of mag- 
nitude less than the numbers we obtain below and shows, 
that even objects billions times less massive than the Sun 
could introduce noticeable variability provided they con- 
tribute to the overall compact object density (however, the 
time scale of this kind of variability will be much shorter). 

The second assumption - which relates 7 and R - is 
of a rather technical nature and does not restrict our scope 
whenever marginal cases of zero shear are not considered. 

As shown in the Appendix B, under these assumptions 
the integral 1321 1 may be then rewritten in the following form: 



;(k,7,_r) = 



I det J\ 



I{K,l) 



ff3(«,7) 



hxR 



1. 



(34) 



where I(k,j) and 173(^,7), denned by 1B22I and <B13ll . are 
computed numerically. The actual values of (73 (k, 7) do not 
exceed ~ 0.1 — 0.2. Therefore we neglect the weak depen- 
dence of e^ on R putting R — 10 -6 . 

The contour lines of constant e M are shown in Figure 
for the range of parameters k, and 7 present in the observed 
clusters. Computations near the lines 7 = |1 — k\ are unre- 
liable and are therefore ignored on the graph. 



2.3 Dependence on source redshift and smooth 
matter contribution 

The actual values of k and 7 depend on both the physical 
surface density at the point in the lensing plane where the 
image is formed and the scaling parameters E and r o given 
by 122I23H which in turn are functions of the distance pa- 
rameter D 1241 and thus - of lens and source redshifts and 
underlying cosmology. If we consider the redshift of the lens 
zi as fixed, the values of convergence and shear would de- 



pend on the source redshift z B . As S oc D _1 and k oc E" 1 
and the same applies to the shear, k and 7 are directly pro- 
portional to D and can be written in a simple form: 



D{zs) 

D(zS) 



(35) 



with k and y being the convergence and shear correspond- 
ing to a given redshift z°. 

For the currently favoured flat cosmological models 

(A = i-n„) 



D(z 8 ) c 
where 

/(*) = 



f{z 3 )-f{ Zl ) 



= 1 






dC 



y/n (t + 1) 3 + 1 - n„ 



(36) 



As f(z) is a monotonic increasing function of its argument 
and z s is clearly greater than zi, the second fraction in 1351 



h(z s 



D(zs) 
D{z°) 



fW/f(*s) 



(37) 



increases from zero at z 3 — zi through unity at z s = z° to 
some limiting value hoc - determined by zi, z° and Q - 
when z s — > 00. This is somewhat different from considering 
convergence and shear as functions of zi , in which case there 
exists an optimal lens redshift which maximizes the lensing 
parameters. In the case of varying z s the further the source 
is the greater k and 7 are. For instance, in the case of Abell 
370 with zi = 0.37, /loo ~ 1.7 for 0, o = 1 and about 1.5 
when £l is only 0.2 (z° = 1) - that is, k and 7 for far away 
sources are not much larger than for sources at redshifts of 
about unity. 

The behaviour of e M with redshift is evidently more 
complex - as k and 7 slide along the line of proportion- 
ality I35H in K-7 plane the variance first increases from zero 
at its bottom left corner but can then, depending on n and 
7o cross one or two 'zero-signal' lines 7 = |1 — k\. Actual 
z s — 1 convergence and shear in Abell 370 and Abell 2218 for 
which we have det ailed density maps feezecourt et alJll999l : 
iKneib et alJll99o) cover approximately the range present on 
figureQtherefore there is no much point in discussing how e M 
changes with redshift any further, specially since the mea- 
surements of the redshift have been performed for many of 
the potential targets for surface brightness variability obser- 
vations. 

However, it is worth noting the general pattern of 
brightness variability behaviour over the area of some of 
the most prominent candidates for this sort of observations 

- gravitational lensed arcs. These objects often consist of 
two or more sections with critical lines between these sec- 
tions and, in the case where there is no smooth matter,the 
variability will be most easily observed in the pixels further 
away from the critical lines on which det J - and variability 

- vanish. 

When compact objects make up only a limited fraction 
of the lensing matter - which is expected to be the case 

- the situation is more interesting. Let x be the compact 
objects share in the total convergence, so that k' — xntot and 
K c — (1 — x)ntot- Then the effective convergence and shear 
are K — XKtot/\l - (1 - x)Ktot\ and 7 = j' /\1 - (1 — x)ntot\- 
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The factor det J in e^ is then: 

2 _ [|1 - (1 — x)ntot\ - xntot] 2 - -y' 2 



det J = (1 — k) — 7 

1 



|l-(l-:r) Ktot | 2 

(l-i^tot) 2 -!' 2 

[1 - (1 - x)Ktot} 2 1 [1 - (1 - 2x)K tot ] 2 - 7 ,; 



(38) 



the latter alternative determined by whether Ktot is less 
(top) or greater (bottom) than (1 — x)" 1 ^ 1. 

Therefore, the outer lines of the zero variability signal 
(those which correspond to 7 = 1 — n and therefore k < 1) 
are not affected by the addition of smooth matter. Other 
zero line positions depend on the value of x and this de- 
pendence represents a potential means to determine this 
value. We will see that only highly magnified pixels show 
variability detectable with present-day observational tech- 
niques - i.e., those lying near the critical curves, and there- 
fore for the effect to be detectable these curves should not 
coincide with the lines of zero variability. Thus, the condi- 
tion Ktot ^ (1 — x)" 1 or, since Ktot can be determined from 
macrolensing modeling: 



x < 1 — 1/Ktot 



(39) 



is in practice necessary to observe the effect. For axially sym- 
metric clusters, the arcs which form on the second (inner) 
critical curves tend to have radial morphology, i.e. their di- 
mensions along the critical curve - and thus the number of 
highly magnified pixels - are small. 

As an example we have computed the maps of the signal 
e M f or two well-studied clusters - Abell 370 and Abell 2218 
(seelKneib et alll99dlBezecourt et al.ll999tlMetcalfe et alJ 
l2003l and references therein) and present them in figures [5] 
and [3] These are given for two values of the source redshift 
z 3 — 2 for both clusters and z s = 0.724 and z s = 0.702 for 
Abell 370 and Abell 2218, respectively. The latter two corre- 
spond to giant gravitational-lensed arcs seen in the clusters 
while the former are given for comparison. As new instru- 
ments - like James Webb Space Telescope (JWST, formerly 
known as NGST) - come into operation they are expected 
to observe many more lensed galaxies behind these clusters 
and z s — 2 maps show how the signal might look for them. 
Each of the maps is given for two values of x with 100 and 20 
per cent of convergence contained in compact objects. These 
values of x are assumed to be constant over the maps. 

In Figure [I] we show the contour lines of the microlens- 
ing signature e M superimposed on the optical image of the 
radial arc R in Abell 370 obtained wi th the Hubble Space 
Telescope (HST) jBezecourt et alll999l). The source redshift 
her e is estimated to b e z s w 1.7 fsee iBezecourt et al.lll999] 
and lSmail et al.lll996l) . The regions between the thick white 
lines correspond to 'zero signal' lines where e M ^ 0.3, while 
the regions between the thick black lines - where present - 
have \k c — 1| <C 1, where the analysis given in this paper is 
not applicable. Dashed white lines show the location of the 
critical curve. 

The figure illustrates how the signal changes with vary- 
ing fraction of compact objects in the overall mass budget. 
Perhaps in contrast to naive expectations, the signal gen- 
erally increases when the density of compact objects is de- 
creased because of the magnification effect by the smooth 
matter distribution. This can be understood on the basis 
of 1341 1: the dependence of the variability extent on the 



source size is rather modest while slight changes in the 
smooth matter convergence k c change the |1 — k c | _1 factor 
in the definition of effective shear and convergence of the 
compact matter IKavser et al. 1986; Paczvhski 1986) signif- 
icantly when re c is about unity which is a common place for 
macrolensed images of distant sources; thus convergence and 
shear on the n — 7 plane of Fig0 can assume high values. 
This somewhat surprising behaviour has also been discussed 
bv ISchechter et aU (2003) while ISchechter fc Wambseansa 
(2002) give a detailed explanation of this effect. More im- 
portant is the change in zero signal lines pattern that can 
be readily probed in observations and can provide, via 13811 . 
an interesting constraint on the compact object contribu- 
tion to the overall convergence determined by modeling of 
the lensing potential. 

For images which do not lie on the critical lines ob- 
servations can still be of interest for the determination of 
x through studying the variability pattern in greater detail 
and comparing it to the predicted one. However, for the 
two clusters investigated in this paper the latter possibility 
remains mostly a theoretical one because of observational 
limitations. 



3 OBSERVATIONAL ASPECTS 

Let us now discuss the prospects for the detection of the 
considered effect. We will consider observations with Hubble 
Space Telescope as a reference point in this section although 
it will be clear that observations of this effect with HST 
in the two clusters under investigation is impractical. Ob- 
servations with more advanced instruments, such as James 
Webb Space Telescope or the proposed 30-metre telescope 
(also known as CELT) 3 could, however, be used to observe 
the microlensing- induced variability. 

The number of photons I detected in a unit time in- 
terval in a pixel from a source of observed luminosity L bs 
(uncorrected for lensing magnification) is determined by the 
luminosity distance to the source Dl(z„), the energy distri- 
bution in its spectrum fx (such that J d\fx = 1) and the 
telescope efficiency rj\ and diameter d: 



1 = 



= L b s d 2 ri e ff\ e ff/16D'i(z s )hc = aL t B 



(40) 



These photons will be accompanied by b background 
photons. For HST Wide Field and Planetary Camera 2 
(WFPC2) and actual luminosities observed in pixels of grav- 
itationally lensed galaxies, the noise is dominated by Pois- 
sonian fluctuations in count numbers. At z s = 1.7 the coef- 
ficient in 140H is a ~ 5 • 10 -6 Lg 1 • hour -1 while the back- 
ground level is about b — 1000 photons in a pixel per hour 
(this value changes by about 1.5 — 2 depending on the source 
heliocentric ecliptic longitude). 

Let us now calculate the time t required to detect the 
fractional change of /3esl in a pixel with given microlensing 
parameters e M , 71 at a signal-to-noise level Q; f3 here deter- 
mines the fraction of pixels deviating from the mean at /3esl 
level via the normal law 1 — $(/3) and we will use /? = 1 for 

3 See http://tmt.ucolick.org/ior details 
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Figure 2. The map of microlensing-induced variabil ity parameter e u over th e area of the cluster Abell 370. Coordinates are given in arc 
seconds, orientation of the images as the same as in lBezecourt et alj J1999I) - the north is to the top, the east is to the left. The maps 
are given for two redshift values z = 0.724 and 2 = 2 and for two values of compact object mass fraction x = 1 and x = 0.2 (for x = 1 
white lines of e M coincide with the macrolensing critical lines, according to 1381 ). The thick black lines at x = 0.2 correspond to regions 
with \k c — 1| -C 1 where the analysis given in this paper is not applicable. 



numerical estimates which is close to the optimal value. The 
signal S is 

S = (3e SL lt (41) 

while the noise N is determined by Poissonian fluctuations 

N = y/zbt + (2 + Peewit, (42) 

the factor of 2 in the latter expression comes from the fact 
that we need to compare two images from different epochs. 
In most cases esl can be neglected for noise estimation. 
Equation 116H shows that signal and noise behave similarity 



which, as can be easily seen, gives the following expression 
for the time required: 



t = ■ 



■ 2(Lj) 1 + b/l 
(L?)a pela 



(43) 



Thus, it is determined mostly by the telescope and geometry 
(through a) and lensing characteristics of the pixel (through 
By. and Jl) while the dependence on photometry is very weak 
as soon as the background value is exceeded by the source 
surface brightness and increases inversly proportional to the 
latter if it is lower than the sky level. In fact, the surface 
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Figure 3. The map of microlensing-induced variability parameter e^ over the area of the cluster Abell 2218. Orientation is the same as 
in figure|5] The maps are given for two redshift values z = 0.702 and 2 = 2 and for two values of compact object mass fraction x = 1 
and x = 0.2(for x = 1 white lines of e M coincide with the macrolensing critical lines, according to 1381 1, The thick black lines at x = 0.2 
correspond to regions with \k c — 1 <C 1 where the analysis given in this paper is not applicable. 



brightness of a typical galaxy at z = is of order 21 m — 22 m 
per sq.arcsec and scales as (1 + z)~' 4_p ' with p depending 
on the spectrum; surface brightness is conserved in grav- 
itational lensing. The sky background outside the Earth's 
atmosphere varies in the range of 22 m — 23 m per sq.arcsec 
and is about half a magnitude higher for the best terrestial 
observatories. Therefore, typical values of the numerator in 
the last fraction of 1431 is of order 1 for nearby galaxies and 
grows rapidly as the redshift exceeds unity. 

Combining values of coefficients in llul and I4UH wc 
find that for observation of the radial arc in Abell 370 with 
HST WFPC2 the value of the first fraction in (Eft is ap- 



proximately 1.2 x 10 4 hours. The value of the variability 
power parameter e^ does not exceed ~ 15 — 20 (in fact, 
e^ changes very slowly with convergence after it exceeds ap- 
proximately 15) while the ratio of sky background to the 
obser ved arc surface brightn ess b/l for most pixels is about 
7-8 llBezecourt et alJll999ft . 

Therefore the effect can be most easily observed in pix- 
els of high magnification JL. This value does not depend on 
the compact-to-smooth convergence ratio, and peaks at the 
critical curve. The variability e M , on the contrary, follows 
the compact matter distribution and wherever some smooth 
matter is present, can preserve high values at the regions of 
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Figure 4. The contour lines of microlensing-induced variability parameter e M superimpo sed on the optical ima ge of the radial gravita- 
tionally lensed arc (z ks 1.7) seen in Abell 370 obtained with the Hubble Space Telescope iBezecourt et aLHl999tl . The area between the 
solid white lines correspond to nearly zero signal with e„ $J 0.3. The dashed white lines represent the approximate location of the critical 
curve, while the area between the thick black lines correspond to \k c — 1| <C 1 where the analysis given in this paper is not applicable. 
Four values of compact to total mass density ratio are assumed: x = 1.0, x = 0.4, x = 0.2 and x = 0.1. 



high magnification. As can be seen from IJ380 this is the case 
when the local convergence value is greater than the inverse 
of the smooth matter share (1 — x)~ , or x ^ 1 — l/i^tou 
otherwise the variability zero lines coincide with the critical 
curves and it is not possible to get both appreciable vari- 
ability and high magnification values. Somewhat ironically, 
the compact objects can only be observed when their mass 
contribution is sufficiently low. 

Magnification values are determined firmly by the 
present-day advanced methods of mass distribution mod- 
eling in lensing clusters which proved to be accurate as 



well as highly and successfully predictive llKneib et al.lll993 : 
Ebbels et al. 1998). The values of convergence for the radial 
arc R in Abell 370 span a range of approximately 1.3 to 1.4 
and therefore the maximum values of the fraction of compact 
matter which produce detectable signal would be approxi- 
mately 23 to 30 per cent - values close to those suggested 
by studies of the Galaxy and its immediate neigh bourhood 
JAlcock et alJl2000l:lLasserre et all200oT:ISadouletll999l) . All 
other gravitationally lensed objects in the cluster either lie 
on the outer critical curves (e.g, the giant gravitationally 
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lensed arc A0) or do not show sufficient magnification val- 
ues. 

For the map of variability present on Figure|l|with x = 
0.2, the values of e 2 ^ on the arc are about four while the 
magnification varies from a few dozen to a few hundreds with 
a handful of pixels where JL exceeds f 4 . Therefore, according 
to eqn. 14311 . for most pixels detecting variability at signal- 
to-noise ratio Q — 5 with HST would require considerable 
integration time of a few hundreds to a few thousands hours. 
However, for a few pixels these exposure times will have 
more reasonable values of order ten hours. 

Certain reservations should be made reflecting the fact 
that these values are dependent on the model and in this 
respect the distribution of the magnifications (or the derived 
required exposure times) is a more robust measure. However, 
one should keep in mind that the critical curve is necessarily 
a set of points with infinite magnification and therefore the 
number of variable pixels is determined by the length of the 
arc along the critical curve (which is rather small for radial 
arcs) and the rate at which magnification falls off the critical 
curve - i.e., graduality in convergence and shear values over 
the image in the critical curve vicinity. 

To estimate the latter value, we can rewrite l|43|l : 



t — to X 



where 



to 



2 2(£i) 1 + b/l 



(44) 



(45) 



which is about 6 x 10 hours for the radial arc R observed 
with HST. Factor 1//Z vanishes at the critical curve, and to 
the first order approximation, t as a function of the coordi- 
nate d orthogonal to the critical curve is 

t = to\V=\d. (46) 

M 

Hence, the width of the strip along the critical line on which 
the required integration time t is less than a given value T 
is simply 



2d: 



2T 



T 



1 



tJVil t |l-«t t||V«t t|+7'|V7'l 



(47) 



For the patch of the critical curve near the radial arc with 
e£ « 4, 1 + b/l « 7 - 8, Ktot -1 = 7- 0.3, |V« to t| ~ 
2 x 10~ 3 pix~\ V-y] ~3x lO^pix -1 and /3 = 1 one gets 
2d « T/600 h for HST. Multiplied by the dimension of the 
arc along the critical curve (about five for the radial arc R 
in Abell 370), these give the required number of pixels. 

It is immediately clear from this estimate that since 
integration time of more than 100 hours is hardly possible, 
only local (on inter-pixel scale) stationary points in 1/71 can 
give detectable signal for HST images - the above mentioned 
strip itself is too narrow. One would need more advanced 
telescopes to observe the effect, such as JWST or CELT, - or 
explore other lensing clusters where strongly lensed objects 
on inner critical curves are seen. 

At visible and near-infrared wavelengths, the sky back- 
ground level expected to be observed with JWST is not 
much different from that with HST and therefore changes 
mostly come from differences in the optics and spectral band 
via the value of a in eqn. 1431 . Using the JMS sensitivity 
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Figure 5. The number of pixels with detectable microlensing 
signature as a function of integration time required for the one 
e M variation to be detectable at five sigma level. 



calculator 4 we estimate that the exposure time required for 
the Near-InfraRed Camera of JWST to detect the signal 
will be about 15 to 20 times as short as those of HST. How- 
ever, even this would require exposure times of several dozen 
hours as 2d w T/30 h in this case. In light of the recent dis- 
cover y of a candidate z = 10 lensed galaxy behind Abell 
1835 JPello et al J 120041) . it does not seem implausible that 
such ultra-deep exposures with JWST will be attempted. 
For the radial arc in Abell 370 that would result in about a 
dozen variable pixels. 

The proposed 30-metre telescope would make the 
prospects more optimistic. For a ground-based telescope, 
the sky background will be a factor of 1.5 — 2 higher, and 
the atmosphere transparency should be taken into account. 
However, increase in collecting area over the HST will be 
enormous, and the net effect will be to reduce t to about 
8 x 10 3 hours. Of further advantage would be the use of 
diffraction-limited mode. Taking into account possible track- 
ing uncertainties, an estimate for the angular resolution of 
~ 0.01 arcsec should be considered conservative. This rep- 
resents a five-fold decrease in the pixel size resulting in the 
ability to get closer to the critical curve. Thus, according 
to 1471 . the width of the strip around the critical curve where 
variability can be detected in T hours integration time would 
make 2d « T/l h . Similarly, the arc length along the critical 
curve will be covered by five times more pixels compared 
to HST. This means dozens of variable pixels and potential 
to observe the pattern of variability change both along and 
across the critical curve in just a few hours long exposure! 

Figure |5] presents the histograms of the number of pix- 
els in the image of the radial arc seen in Abell 370 which 
are expected to show variability detectable at the 5a level 
as a function of integration time for JWST and CELT. The 
actual values of the arc brightness are used to determine t 
from H43t for every pixel which are then binned logarithmi- 
cally in 0.5 dex wide bins. Two values of the compact matter 
fraction x = 0.2 and x — 0.1 are assumed. For JWST, the 

4 See http://www. stsci. edu/jwst/science/jms/jms_flux_form. html 
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pixel size and background level are taken to be equal to those 
for HST WFPC2. For CELT, the background level is half a 
magnitude brighter while for the pixel size a value of 0.01 
arcsec is assumed. 

One should bear in mind that the pixels mentioned 
above are variable at a detectable level, but they will only 
spend about one-third of the time in this 'varied' state. This 
fraction of time can be controlled by the parameter j3 but 
one can see - combining 14311 . the normal distribution and 
the nearly linear histogram shape in the region of interest - 
that the value j3 = 1 is close to optimal. 

The time scale of variations depends on the microlens 
masses and motion and is of order months to years for Solar 
masses mov ing with a velocity of a few hundred kilometres 
per second IlLewis e t al. 2000). Observation epochs should 
therefore be separated by a similar time interval. 



4 DISCUSSION 

We have shown in the previous section that with a few ex- 
posures on JWST or CELT we expect to detect shimmering 
of pixels in the image of the radial gravitationally lensed 
arc in Abell 370 due to microlensing by compact objects in 
the cluster. To answer the question of how such detections 
should be interpreted let us now sketch a portrait of a typi- 
cal event. This would also give us an insight into what sort 
of contaminants could mimic variability due to gravitational 
microlensing. 

The surface brightness of the arc is about seven to eight 
times lower than the sky background level outside the Earth 
atmosphere. In a 30 hours-long exposure on JWST about 60 
thousand 'signal' photons will be detected in a typical pixel 
and around 500 thousand background photons will accom- 
pany them. Hence the noise, according to I42H . is approxi- 
mately (neglecting ^ssl) 1060 photons and five times that 
is 5300 photons or about 9 percent of the original flux. As 
Sft in this region is about 2 (see Fig0Jl!), the pixel, according 
to 11611 contains light of approximately 4.7 x 10 4 stars with 
an intrinsic luminosity of some 1.4 x 10 4 Lq. Given a typical 
observed luminosity of 2 x 10 s Lq the magnification needed 
is ~ 10 4 , in accordance with the estimate of the previous 
section. 

Nine per cent variability corresponds to 0.09 x 2 x 
10 8 //II/q = 1.8 x 10 7 /TIL© and the average value of magnifi- 
cation in the map is about a hundred. In the case of the aver- 
age pixel, only supernova and brightest peak nova eruptions 
can give the true increase of ~ 10 5 Lq. The most significant 
contaminant to the average pixel is nova eruptions. We can 
calculate the ex pected number of n ova explosions in a way 
similar to that of B altz et alj J2003T) - for a galaxy similar to 
the Milky Way the rate of nova explosions is expected to be 



are not a problem for our study. Clearly, supernova expo- 
sitions in the source galaxy are even less of a problem. An 
additional source of contamination is supernovae in back- 



ground galaxies but, with a rate of 10 



10" 2 SNe per 



10" 



■10" 



eruptions per star per year. Eruption durations 



- by which we mean the period of time novae stay above the 
level of interest ~ 10 5 - is only a fraction ~ 0.01 — 0.1 of 
a year, long compared to the integration time while short 
compared to the interval between exposures. Therefore we 
expect about 10 -10 erupted novae per star at any given ex- 
posure. The radial arc spans around three hundred pixels 
and therefore the expected number of stars in it is about 
10 - i.e., this is a galaxy of rather modest size. Thus, even 
not taking into account the multiple nature of the arc, novae 



galaxy between two exposures (see Sarajedini, Gilliland and 
Phillips 2000), they are not important. 

However, with pixels magnified by a factor of ~ 10 4 
which is needed to observe the microlensing variability, phys- 
ical luminosity changes go down to about two thousand Solar 
luminosities or less, and this is about the amplitude of the 
brightest Mira variables in red bands. Contamination due to 
Miras (as well as other variable stars) in microlensing stud- 
ies is usually removed by considering observations in differ- 
ent spectral bands. The three key signatures of microlens- 
ing origin of the variability are achroma ticity, uniquenes s 
and a symmetric form of the variations JPaczvriskil Il986) . 
None of them is valid in the case considered in this paper. 
Achromaticity does not work for pixel lensing although in 
the case of low optic al depth some constraints can still be 
applied llGouldlll996l) . When the optical depth is high every 
star in the pixel at any given moment is subject to strong 
microlensing and this fact does not allow to use achromatic- 
ity constraint. For the same reason the uniqueness of the 
microlensing event does not work any more. With regards 
to the symmetry, we do not observe individual light curves 
in this case and therefore cannot use this constraint at all. 

However, although we cannot use achromaticity for in- 
dividual pairs of measurements, this property is still valid in 
a statistical sense. Namely, the variability extent seen in dif- 
ferent bands is, according to 11611 . proportional to {L 2 ) /{Li) 2 
- the value which does not change much from one band to 
another, is closely related to the magnitude of the surface 
brightness fluctuations and can be determined in observa- 
tions of nearby galaxies which are definitely not lensed. In 
contrast, (absolute) variability amplitude of variable stars is 
strongly dependent on the spectral band and, for instance, in 
Miras the change ranges from thousands Solar luminosities 
in K band to hundreds and even tens of Solar luminosities 
in bluer bands. Other variables are too faint to affect the 
fluxes of pixels containing thousands of stars. 

Contamination due to variable stars is a more serious 
issue for observations with CELT. An analysis similar to the 
one given above, shows that typical intrinsic luminosities of 
pixels with variability detectable in a one hour long expo- 
sure amount to just around 500 solar luminosities with typ- 
ical variability extent of around 70 per cent or around 350 
solar luminosities. This is a range at which various variable 
stars may contribute to the observed variability. The only 
way to distinguish it from the variability due to gravitational 
microlensing by compact dark matter in the cluster seems 
to be to use the behaviour of variability from one pixel to 
another and across spectral bands to see whether it is consis- 
tent with physical variability or gravitational microlensing 
hypotheses. 

One of the remaining problems is how to tell the dif- 
ference between no compact matter and too much compact 
matter in the case of a null signal detection. More work 
needs to be done on this question and perhaps other ef- 
fects should be considered to answer it. However, the effect 
considered provides us with a lot of information on the mi- 
crolensing population. The Gaussian approximation seems 
to be a simple framework for characterizing pixel microlens- 
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ing in galaxy clusters and although the implications of future 
observations on microlensing population are not straightfor- 
ward, they can provide strong constraints on the otherwise 
unaccessible properties of this population. 
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APPENDIX A: DERIVATION OF 
FORMULA (30) 

The deriv ation pres e nted here follows very closely the lines 
of that bv iNein dorf (2003J). We do not include all the steps 
of this derivation which can be found in the original work. 

Starting with the expression for the magnification fac- 
tor 1271 one has the following result for the observed flux of 
a source with surface profile /'(£ — £ ) = I /(ttR 2 ) x 0(_R — 
C - C°l) placed at £ : 



/(Co) 



1 — K c 



d 2 Cd 2 zJ'(C-Co)5(C-^z + S(z))(Al) 



Its average value over S(z) is well-known to be independent 
of microlensi ng population mass distribution a nd source pro- 
file (see, e.g. lSchneider. Ehlers fc Falcolll99l Chapter 11) 



(I) 



la 



|[l- Kc H(l- K )2- 7 2] 
Io 



(A2) 



|l-K c | 2 |det J\ 
where J = J — rel. 



The value of I 2 (Co) is calculated in a similar manner: 
/ 2 (Co) = h * n /d 2 Cid 2 C 2 /'(Ci - Co)/'(C2 - Co)x (A3) 

1 Kc\ J 
d 2 Zld 2 Z2<5(Cl - JoZl + Sl)5(C2 - JoZ2 + S 2 ) 



5 Please note, however, that the asymptotes for the involved func- 
tions found here differ slightly from those found in JNeindorfl 
2003) due to some numerical error and a typo in the latter work. 
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tt/2 



where Si = S(zi), Si = S(zi). Introducing the joint proba- 
bility function density of Si and S2 y(Si, S2, zi, Z2) we can 
calculate the average value of I : 



(n 



d 2 Cid 2 C2d 2 zid 2 z 2 d 2 Sid 2 S 2 (A4) 



1-K C | 4 

S 2 xC 2 xTZ 4 

5(Cl - J0Z1 + Sl)5«2 - JoZ2 + S 2 ) X 

^(Si,S 2 ,Zi,Z2)/'(Cl)/'(C2) 

Changing to the Fourier domain and making use of S - func- 
tions we get: 



<I 2 > = 



(2tt) 4 |1-k c 



d Tld T2d Zld Z2 



(A5) 



TZ 2 xC 2 



/(Tl)/(T2)Q2(Tl,Zl,T2,Z2)e- l(TlJ ° Z1+X2J "° Z2) . 

where 

Q2(Ti,Zl,T2,Z2) = 

J d 2 Sid 2 S 2¥ ,(Si , S 2 , zi , z 2 )e^ lSl+ ^ S2) 

R 1 



(A6) 



is the Fourier transform, or characteristic function, of 
y(Si, S2, zi, Z2) and I(t) is the one of the source profile, 
which is the following for a uniformly radiating disc of ra- 
dius R emitting a total flux I : 



I{r) = ^MRr) 



(A7) 



We introduce central and relative coordinates t\ — T — 1/2, 
t 2 = T+t/2, Zi = R c r/2, z 2 = R c +r/2, assume that lens 
positions are not correlated and take the limit of an infinite 
lens plane to obtain the following expression (J = J ): 

1 1 



(n = 



(47r) 2 |l- Kc | 4 |detJ| 



(A8) 



d 2 td 2 r7(-i)j(|)Q( t) r) e -i^ Tt > : 



where the 5-functions in R c and T have been utilized and 
for 

Q(t,r) S Q a (-i > -|,i > J) 

the following expression was obtained (|Neiridorj[2003) 

Q(t,r)=e-* r2a(s) . (A9) 

Here s — - and the angle x between s and positive Ox ray 
equals the angle between t and r. Function a(s,x) is the 
mass average of a(ms, \) 



a(s,x)= / dmcf>(m)a(ms,x), 
the latter given by 



1 — e lsx 
a(s,x)= I dx — —^ — f(x,x) 



Finally, the function f(x, \) is defined as 



(A10) 



(All) 



cos (f>dcj> 



-tt/2 



y(x — 2 cos x) 2 c °s 2 4> + ( x sin + 2 sin x cos < 



and can be expressed analytically in terms of complete el- 
liptic integrals: 

1 — x cos x 



f(x, X ) 



u + 1 I 



1 



K{v)+ 



(A12) 



(1 — a;cosx)(l + u )' 2 

2~U? 



[E(v)-K{v)]} 



with u — \/x 2 — 2x cos x + 1 an d V — 4u/(u + l) 2 . 

The function a(s, x) is therefore easily computed nu- 
merically. However, its overall behaviour is easily guessed 
from the following two analytic asymptotes: 



Oi(s,x) "^ irs-i- cosx (A13) 

valid to the accuracy of ~ 1 per cent at s ^ 3 — 5 and 

a(s,x) S -^ Js 2 (i+21n2-7 + cos 2 x-lns) (A14) 

— ins cosx 
fa — s (l.3 + cos x ~ Ins) — in cosx 

n 

where 7 = lim (^ 1/fc — Inn) (a 0.577216 is Euler's con- 

n -"°° fe=i 
stant. 

Substituting the expression 1A7I for I(t) into !A8t . ex- 
panding 

--(J D t)r = --rt[cosx + 7cos(x + 2a r )] 

and using the change of variables t = rs we can immediately 
integrate over the angular component a r of r to get 

2/„ 2 



{I) = — 

7r|l-K c | 4 |detJ|i? 2 

00 00 2tt 

dr/d S /dx^i 2 (^) 



(A15) 



j ( jLJL \ e -T[ n< " - x : - "- \ 



Thus, using expression <A21 for the average value of ob- 
served flux we obtain the integral 13UII . 



APPENDIX B: EVALUATION OF 
INTEGRAL (32) 

To evaluate the integral in (I32H we consider the asymp- 
totes of B(a,x) which follow directly from the asymp- 
totes 1A14I lA13ll of the function a(a, x) introduced in the 
Appendix A: 

t-., v K ,, „ 2 (m 2 lnm) , . , T , H . 

B(cr, x) « ^rn cff a(1.3 + cos 2 X - V , 2 . ' - In a) (Bl) 
2 [m ) 

+i(l — «:)cosx, J«l 

B(a,x)~B( X ) = K + icosx, <r » 1 (B2) 

with m cft = (m 2 )/{m) jRefsdal fc Stabelllll99lt) . 

At 7 ^ 1CF 2 it is convenient to split the integration 
over s into the following five regions: 

(i) 00^ crii? 2 /2, 



(ii) criir/2 sC s < a 2 R 2 /2, 



a x < 1 
cr 2 > 1 
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(iii) a 2 R 2 /2 < a < (J37 2 , 
(iv) 0-37 2 < s < (T4, 
(v) (74 < S 



(73 < 1 
CT 4 > 1 



Due to the assumptions on 7 and R made, in regions (i - iii) 
the convergence of the integral in p is provided by Jo(yp) 
and e~ Bp while J 2 (^/sp) w sp/4 holds well for all values of 
p where the integrand is any significant. Therefore 



dpJ 2 ( v / sp)Jo(7P)e 



-pB 



dppj (-yp)e 



B 



(B3) 



4(B 2 +7 2 ) 3 / 2 ' 

In addition, in regions (i) and (iii) approximations IBlll 
and ll I i> 12 D can be used, respectively. The latter approximation 
is even better in regions (iv) and (v) but (113311 clearly fails 
there. Let us consider the five regions in turn: 

1. OSC S < CTl -^ 

Returning to the variable a — 2s/ R 2 we can write down 
this part of the integral as follows: 



B(cr,x)dx 



[B 2 ((7, X )+7 2 ] 3/2 





We can rewrite B of IBlll in the following form: 

B{cr,x) = \1- K\(x±icosx), 

where 

_ K m e « , 2 (m 2 lnm) 

X = — — r(1.3 + cos x — — ; — 57 — — maja 



2 11 — as: 



<m 2 ) 



(B4) 



(B5) 



(B6) 



and at sufficiently small a when the logarithm term domi- 
nates x is nearly independent of X- Then 



h 



(1(7 



(7 \\-k\- 



5i x(a), 



7 



1-K 



where 
gi(x,5) = 2 



(x + i cosx)dx 

[5 2 + (x + icos X ) 2 f /2 



(B7) 



(B8) 



We have found numerically that the following approximation 
holds for gi(x,5) with an accuracy of ^ 1 per cent: 



1 x\ 5i(l 1 -^l.< 5 ) 

gi(x,5) ~ — - — j — x, 



X^Xmax = \l-5\ (B9) 



Xmax stands for the point where g\ reaches its maximum at 
a given 8. This maximum value is approximately |1 — S^ 1 for 
|1 — S\ ^ 10~ 2 and should be calculated directly otherwise. 
At |1 — 5\ > <7i lnai we can therefore perform an inte- 
gration over a in 1B7> to get 



h 



1 /tm. 



<n 



1- 



421- 



1.8 



1 — k\ — 7 j 



(BIO) 



{m In m) 

(m?) 



In — 1 (7i 

(71 



When (7i does not obey the condition just formulated, 



we just replace it with a[ < o\ such that a[ In a[ is smaller 
than |1 — 7/|l — k|| and move the rest of the calculations to 
region (ii) where it is performed numerically. 

The problem arises with |1 — 8\ — > but this corre- 
sponds to the case of diverging average amplification and 
the microlensing-induced variabil ity is expected to drop log - 
arithmically to zero in this case llDeguchi fc WatsorJll987f) . 
Formally, it does happen in our calculations - det J in front 
of the integral 132H make it go to zero linearly in all regions 
but (i) where the ratio of the determinant and denomina- 
tor of IB 1CH tends to a non-zero limit, while gi&[ lnaj 
provides the behaviour expected. However, numerical cal- 
culations become unreliable in this case and so we do not 
calculate e 2 for (k, 7) closer than approximately 0.03 to the 
7 = 1 — re I lines in re-7 plane. 

Another apparent problem with (IBlOl l seems to be 
present when re — » 1, but this turns out to be a slight techni- 
cal issue with no real computational consequences and there- 
fore may be called 'a removable discontinuity'. 

2. cjiR 2 /2 < s < a 2 R 2 /2 

This region is the easiest to compute: 



0"2 



where 



52(0") 



B(a, X )d X 



[B 2 ((7, X )+7 2 ] 



213/2 



(Bll) 



(B12) 



and B(a,x) is computed numerically by interpolation of 
a(a, x) which is computed in advance with good accuracy. 
The computations of I2 provide no problem from either 
conceptual or computational point of view but it transpires 
that they are the most time consuming part of the proce- 
dure. 

3. a 2 R 2 /2 ^ss^ (7 3 7 2 

At s ^ a 2 R 2 /2 the function B(2s/R 2 ,x) - re + icosx 
to high accuracy and does not depend on s. Therefore the 
integral over x : 



53 («, 7) 



(re + icosx)dx 
[(re + icos X ) 2 +7 2 ] 3/2 



(B13) 



(which equals g\ introduced above for x = re and (5 = 7) 
turns out to be a common factor and the integral with re- 
spect to s is elementary: 



1 2(737 

-fl 3 («,7)ln -j 

4 <7iRr 



(B14) 



4 & 5. s> a 3 ~f 2 

In regions (iv) and (v) the approximation 1B2II is still 
valid, therefore 



oc 00 

/j - /.-, - / -j / dpJ 2 ( v ^5)Jo(7p)x 



(B15) 



<?3~f 



II 
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dxe 



-p(ft+icos x) 



= 2tt / —g 4 (s,K,-y), 



where 



g 4 (s,K,-y)= / dpJl( y /sp)Jo('yp)Jo(p)e Kp . 



(B16) 



The integrand in g 4 decreases exponentially and therefore it 
is sufficient to perform a numerical integration up to some 
Pmax/n. The absolute value of the residual can be easily 
estimated 



fl = 



<Wi (y/sp) Jo(7P) Jo(p)e~ 



Ape Kp — —e p " 



(B17) 



The value p m ax ~ 30 turned out to be suitable for all our 
purposes. 

For completeness we can write down the integration in 
region (iv), which is done numerically: 



with y/Tp/2 - in the whole region of actual dependence of 
B(2s/R 2 ,x) on s that makes the result virtually indepen- 
dent of R allowing it to appear in region (iii) only. 

Considering the constants <n used in the actual cal- 
culations we found that the following values provide the 
best compromise between the accuracy of the computations 
and the time needed to perform them: <j\ = 0.03, 02 = 4, 
(T3 = 0.03, (J4 = 30. This corresponds to the case when 
<j>{m) — 8{m — 1) which is the case actually considered. We 
have checked that changes in these constants do not affect 
the results. 



h = 2tt 



As 

— g 4 (a,K,7). 



(B18) 



The integration in region (v) is accomplished by considering 
the asymptotic behaviour of 54 at large s. The integral 1B16II 
effectively splits into two - with p ^ Po/s and p ^ po/s with 
p ~ 1. Invoking the asymptotics of Bessel functions one 
can see that the second of these integrals is proportional 
to l/\/s and represents a leading term when s —* 00. Using 
the asymptotic formula Ji(x) ~ y/2/(nx) cos(ivx/2 + a) and 
noting that the cos-term oscillates rapidly for x>lwe find 
the following limiting value for g 4 : 

g 4 (s, k, 7) -> — t=35(«, 7) , 



where 



TVy/S 



dp 



fls(«,7)= / —=Mp)Mlp)e~~ 




(B19) 



This integration is again done only up to some p m ax/K- 
Thus we obtain the last portion needed to compute e 2 - 



h = 27T 



ds 4 g 5 («,7) 

— g4(a,tt,7)«- 3/2 



(B20) 



Putting all the pieces together we may now write the 
result: 



1 , ^ I dot ./ 



where 



I^)- 9 -^lnR 



Tl \ - T I T I T I r I 93{^l) , 2CT37 2 

7(k, 7) = ii + i 2 + J4 + is H In 



(B21) 



(B22) 



We should note here that it is exactly the behaviour 
of Ji(y/sp) - namely the possibility to approximate it 



